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摘 要 : 基于 祁连山 西 段 老虎 沟 流域 断面 的 径流 数据 与 大 


营 气 象 站 的 气象 资料 ,对 冰川 区 径流 与 


气象 要 素 之 间 的 相关 性 进行 了 分 析 , 并 建立 了 多 元 指数 非 线性 回归 方程 对 径流 进行 重建 以 弥补 缺 


失 径 流 , 此 外 ,对 冰川 区 径流 的 年 际 、 季 节 、 日 尺度 的 变化 特征 进行 了 分 析 。 结 果 表 明 :(1) 冰川 


Ba 


径流 量 与 气温 的 相关 性 最 高 0.86, 其 次 是 水 汽 压 (0.81)、 相 对 湿度 (0.46)、 降 水 量 (0.27) ,径流 受气 


温 影响 最 大 。(2) 21 世 纪 观 测 资料 显示 日 平均 径流 


量 为 2.10 ms"', 较 20 世 纪 50 年 代 末 的 1.65 mts 有 
所 增加 ,主要 因 消 融 季 气温 上 升 0.75 所 致 , 强 消融 期 (7、8 月 ) 径 流量 的 年 际 变化 较 大 ,消融 期 初 
(5、6 月 ) 和 和 末 (9 月 ) 年 际 变化 较 小 。 消 融 季 5 一 9 月 产 流 占 比分 别 为 5.3%、16.1%、37.3%、35.1%、 
6.2%。(3) 多 元 指数 非 线性 回归 方程 可 以 较 好 地 模拟 日 径流 量 (纳什 效率 系数 平均 为 0.70) ,补充 缺 
失 径 流 后 ,对 于 径流 的 日 变化 ,消融 期 初 和 末日 变化 较 小 ,而 强 消融 期 径流 的 日 变化 较 大 。 对 于 径 
流 的 时 滞 效 应 ,老虎 沟 流域 消融 季 各 月 径流 哇 现 “ 谷 - 峰 ” 的 日 变化 特征 ,消融 期 初 和 未 气温 达到 最 


大 值 与 径流 达到 最 大 值 的 时 间 间 距 较 大 ,而 强 消融 期 较 小 ,在 消融 期 初 的 6 月 两 者 相差 最 大 为 


3h。 明 晰 冰川 融 水 变化 特征 对 干旱 区 水 资源 管理 ,生态 保护 及 社会 经 济 可 持续 发 


关 键 词 : 径流 重建 ; 气候 变化 ; 径流 变化 ; 老虎 沟 流 域 
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冰川 作为 “固体 水 库 ” 是 河流 补给 与 农业 灌溉 
的 重要 水 资源 来 源 " ,同时 是 反应 气候 变化 的 天 然 
记录 器 和 预警 器 中 。 冰 川 及 其 融 水 是 中 国 西北 干旱 
区 重要 的 水 资源 ,对 维系 本 地 区 生态 平衡 与 社会 经 
济 可 持续 发 展 具有 重要 的 意义 ”"。 全 球 气候 变 暧 ， 
已 使 全 球 大 多 数 冰 川 出 现 萎缩 。 祁 连 山西 段 自 小 
冰期 至 1956 年 冰川 面积 减 小 幅度 为 16.9% ,冰川 储 
量 减 少 14.1%,1956 一 1990 年 冰川 仍 以 退缩 为 主 ,此 
时 有 段 冰 川 面积 和 储量 减少 量 占 1956 年 的 10.3% 和 
9.396^', XE ET 12 号 冰川 是 祁连山 区 最 大 的 山谷 
冰川 , 自 1957 年 以 来 呈 退 缩 状 态 ,1960 一 2015 年 末端 


收 稿 日 期 : 2022-05-02; 修订 日 期 : 2022-08-15 


具有 重要 


后 退 了 402.96m, 长 度 从 10.1 km 减少 到 9.7 km *' , YK 
川 变化 势必 对 干旱 区 水 资源 造成 重要 影响 。 随 着 
冰川 变化 ,冰川 融 水 对 水 资源 的 补给 也 发 生 了 变 
化 ,这 一 变化 呈现 “ 先 增 后 减 ” 的 趋势 ,如 在 叶 尔 羌 
河流 域 模 拟 发 现 自 1968 一 2017 年 以 来 ,冰川 径流 总 
体 呈 增加 的 趋势 , 且 占 总 径流 的 比重 为 51.1%D)， 
“ 先 增 后 减 ” 的 拐点 还 未 出 现 。 冰 川 径流 “拐点 ”出 
现时 间 与 流域 冰川 大 小 和 多 少 有 关 ,冰川 覆盖 率 
低 、 以 小 冰川 为 主 的 流域 ,拐点 已 出 现 ,如 石 羊 河流 
域 中 ,天 山北 坡 玛 纳 斯 河流 域 和 呼 图 壁 河流 域 ,以 及 
青藏 高 原 的 经 江 源 .黄河 源 和 澜沧江 源 “ ;以 较 大 
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冰川 为 主 的 流域 在 未 来 10~20 a 会 出 现 冰川 融 水 拐 
点 ,如 天 山南 坡 的 库 车 河和 木 札 特 河 .祁连山 焉 勒 
河 及 长 江 源 等 %*。 此 外 ,冰川 融 水 也 容易 造成 自 
PRE ,如 冰 湖 演 决 洪水 与 冰川 泥石流 。 因 此 ,了 
解 冰川 区 水 文 变化 特征 ,量化 当前 和 未 来 干旱 区 冰 
川 融 水 的 贡献 对 水 资源 管理 和 灾害 防治 有 重要 的 

冰川 融 水 、 季 节 积 雪 融 水 、 降 雨 和 地 下 水 是 冰 
川 区 径流 的 主要 来 源 "。 然 而 , 产 流 除了 受气 象 因 
素 的 影响 外 ,也 会 受到 下 垫 面 ` 冰 下 排水 系统 .冰川 
表面 河 网 及 冻 土 活动 层 的 影响 。 由 于 观测 条 件 的 
局 限 性 ,对 产 汇流 的 物理 过 程 认 识 是 有 限 的 。 到 目 
前 为 止 ,基于 度 日 和 能 量 平 衡 的 水 文 模型 可 以 模拟 
和 预测 冰川 区 径流 变化 ,如 分 布 式 能 量 平 衡 模 型 分 
布 式 增 强度 日 模型 和 降水 -径流 - 蒸 散 发 模型 。 
然而 ,在 高 寒山 区 由 于 观测 条 件 的 限制 ,缺乏 较为 
详细 的 水 文 、 冰 川 和 气象 数据 ,大 多 数 水 文 模型 的 
实用 性 受到 限制 “。 在 实际 应 用 中 ,基于 气象 要 素 
与 径流 之 间 的 关系 进行 评估 与 预测 径流 是 非常 重 
要 的 。 时 间 序 列 的 相关 性 分 析 是 确定 径流 与 气象 
要 素 之 间 关 系 的 重要 方法 。 相 关 人 研究 在 冰川 区 已 
经 建立 了 径流 与 气象 要 素 之 间 的 回归 关系 并 总 结 
了 利用 水 文 气象 数据 发 展 的 各 种 回归 模型 “” ,这 
些 结果 表明 利用 气象 变量 可 以 成 功 模 拟 冰川 区 径 
流 。 老 虎 沟 流域 位 于 祁连山 西 段 高 海拔 山区 ,自然 
条 件 恶 劣 ,径流 观测 仪器 (水 位 计 ) 在 夏季 易 受 到 河 
道 砾石 撞击 从 而 导致 部 分 径流 数据 缺失 。 此 外 , 老 
虎 沟 流域 的 相关 人 研究 仅 集 中 在 某 一 年 份 的 径流 变化 
人 研究 且 仪 考虑 气温 、 降 水 恢复 缺失 数据 。 

因此 ,在 高 寒山 区 ,建立 径流 与 各 气象 要 素 之 
间 的 回归 模型 对 径流 时 间 序 列 的 延长 ,缺失 数据 的 
弥补 及 不 同时 间 序 列 径流 的 分 析 有 至 关 重 要 的 意 
义 。 此 外 ,以 往 的 研究 侧重 冰川 区 径流 的 年 际 变 
化 ,针对 径流 的 季 闻 日 变化 的 综合 特征 缺少 系统 
分 析 。 本 研究 基于 架设 在 老虎 沟 12 号 冰川 末端 和 
大 本 营 之 间断 面 的 8a 水 文 数据 和 大 本 营 气 象 站 的 
气象 数据 ,采用 了 统计 的 方法 分 析 了 径流 与 气象 要 
素 之 间 的 关系 ,建立 了 多 元 指数 非 线 性 回归 模型 来 
重建 径流 ,以 弥补 径流 缺失 部 分 。 其 次 ,对 径流 的 
年 际 ` 季 节日 变化 以 及 径流 的 时 滞 效 应 进行 了 系 
统 性 综合 分 析 。 


1 研究 区 概况 


老虎 沟 流 域 位 于 青藏 高 原 东 北 祁连山 西 段 的 
mI EWE , WEEE EE N 96°31'~96°33' E 39°25’ ~ 
39°30'N 之 间 ( 图 1)。 流 域 全 年 受 北半球 中 纬度 的 
西风 带 控制 ,常年 低温 ,降水 量 丰 富 ,流域 水 文 断 面 
设 定 在 大 本 营 和 12 号 冰川 末端 之 间 ,海拔 为 4192 
m, 水 文 断 面 控 制 4~14 号 冰川 ,冰川 编号 为 
5Y448D0004- 6Y448D0014 " ,为 大 陆 性 冰川 。 根 
据 第 一 次 冰川 编目 ,老虎 沟 流 域 冰川 面积 为 24.49 
kmz09。 现 有 断面 为 2015 年 开始 投入 使 用 ,直至 现 
在 ,流域 控制 面积 约 37.59 km ,依据 第 一 次 冰川 编 
目 ,流域 冰川 覆盖 率 为 65.15% ; 2012—2014 年 断面 
位 于 现在 断面 的 东 侧 约 5-10 m, 流 域 控制 面积 约 
37.58 km’, 流域 冰川 覆盖 率 为 65.17% ; 2009—2011 
年 断面 位 于 现 有 断面 的 下 游 约 20 m ,流域 控制 面积 
约 37.67 km' , 流域 冰川 覆盖 率 为 65.01%。 水 文 断 面 
变化 前 后 ,其 控制 的 流域 面积 变化 不 大 ,流域 冰川 
覆盖 率 变化 亦 不 大 , 约 为 65%, 故 水 文 断面 处 观测 
的 径流 可 以 反映 冰川 径流 变化 。 其 中 ,老虎 沟 12 号 
冰川 是 祁连山 最 大 的 山谷 型 冰川 ,长 度 为 9.7 km, , 面 
积 为 20.37 km ,该 冰川 的 流向 为 西北 一 东南 ,海拔 
范围 为 4250~5483 m。 近 年 来 随 着 全 球 气候 变 暖 ， 
冰川 已 经 经 历 了 大 量 的 物质 损失 与 退缩 ,在 过 去 
的 50a 间 ,12 号 冰川 退缩 了 403 m, 面 积 减少 约 
1.54 km“ 。 对 于 平衡 线 而 言 , 平 衔 线 的 平均 海拔 
为 4380 m" ,但 它 在 2010—2012 年 已 经 上 升 到 
5015 m?" , 


2 数据 与 方法 


2.1 气象 和 水 文 监测 

老虎 沟 流 域 在 海拔 4180 m ,4550 m 和 5040 m 分 
别 架设 了 自动 气象 站 ,本 研究 选取 海拔 4180 m KAS 
营 处 的 气象 站 进行 数据 分 析 ( 图 1)。 自 动 气 象 站 使 
用 Campbell 数据 采集 器 CR1000 和 型 号 为 HMP45C 
的 气温 传感器 (精度 为 +0.2 9C) ,降水 记录 使 用 
T200B 自动 十 雪 量 计 , 精 度 为 0+0.1 mm, 每 10s 采 集 
一 次 数据 ,30 min 输 出 一 次 平均 值 ,降水 采用 叶 柏 生 
等 中 提出 的 中 国 降水 误差 观测 及 修正 方法 进行 修 
正 。1959 一 1960 年 气温 、 降 水 量 和 径流 量 数据 均 为 
历史 观测 资料 。 水 文 监测 主要 是 对 流速 和 水 位 的 
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图 例 

K 大 本 营 气 象 站 高 程 /m 
@ 老虎 沟 水 文 断面 per 5471 
E 老虎 沟 流 域 冰川 - 4191 


图 1 研究 区 概况 
Fig. 1 Asketch of the study area 


监测 ,流速 采用 重庆 华 正 水 文 仪器 厂 生 产 的 LS25-1 
型 旋 浆 式 流 速 仪 进行 测量 ,水 位 采用 HOBO 压力 水 
位 计 进 行 观测 ,同时 在 断面 侧面 用 水 尺 进行 验证 ， 
从 而 获得 水 位 -流量 关系 曲线 计算 消融 季 的 流量 。 本 
人 研究 采取 的 主要 水 文 数据 为 :1959 年 5 月 1 日 一 9 月 
20 日 .1960 年 5 月 1 日 一 9 月 3 日 .2009 年 5 月 1 日 一 
8 月 13 日 .2012 年 5 月 1 日 一 7 月 6 日 和 7 月 24 日 一 9 
月 30 日 .2013 年 5 月 1 日 一 9 月 30 日 .2015 年 5 月 11 
日 一 9 月 12 日 、2016 年 5 月 11 H—8 H 31 H 、2018 
年 6 月 12 日 一 8 月 27 日 。 
2.2 研究 方法 

采用 自 相 关系 数 与 互相 关系 数 来 分 析 径 流 自 
身 及 其 与 气象 要 素 的 相关 性 , 自 相关 与 互相 关系 数 


计算 公式 分 别 为 : 


1 
=~ 


i (1) 
六 zx- 
» (x. - XJ(Y.., - Y) 
r,= = z (2) 


式 中 : R, 为 自 相 关系 数 ;n 为 时 间 序 列 数据 的 个 数 ; 
kA ae; X; 为 时 间 序列 日 数据 ; X 为 时 间 序 列 日 
数据 的 平均 值 ;X;,; 为 时 清 为 大 的 时 间 序 列 日 数据 ; 
7 为 互相 关系 数 ;下 为 时 小 为 上 的 时 间 序 列 日 数据 ; 
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7 为 时 间 序 列 日 数据 ;了 为 时 间 序 列 数据 平均 值 。 
采用 变异 系数 来 量化 时 间 序 列 数 据 的 离散 程 
度 ,计算 公式 如 下 : 


(3) 


式 中 : C, 为 变异 系数 ;n 为 数据 个 数 。 

采用 纳什 效率 系数 和 相对 误差 来 评估 模拟 径 
流量 与 观测 径流 量 之 间 的 拟 合 效果 ,计算 公式 
如 下 : 


n 


È Cas E Qua 
E,21- 53 (4) 
YQ. E D 
i=1 
(3o... ~ > 
E,=100x <= 二 (5) 


0. 
式 中 ; E 为 纳什 效率 系数 ; BAAR, 0 ， 为 
观测 值 ; 0 为 模拟 值 ; 0 为 观测 值 的 平均 值 。 


(a) 气温 、 降 水 量 年 际 变化 
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结果 与 分 析 


气温 、 降 水 变化 特征 分 析 
气象 观测 中 气温 和 降水 的 观测 对 水 文 研究 非 
常 重要 ,它们 是 水 文 模型 得 以 发 展 的 基础 。 在 高 海拔 
地 区 ,由 于 地 形 复杂 ,环境 恶劣 ,气象 观测 比较 少 人 3。 
然而 ,高 海拔 地 区 作为 江河 的 发 源 地 ,降水 对 径流 
贡献 较 大 ,冰川 区 冰川 融 水 对 气温 也 比较 敏感 。 此 
外 ,气温 与 降水 的 时 空 分 布 不 同 于 平原 地 区 。 
此 ,加 强 高 海拔 地 区 的 气象 观测 对 分 析 气 候 变 化 背 
景 下 山区 径流 量 的 变化 至 关 重 要 。 老 虎 沟 大 本 营 
气象 站 观测 表明 ,2010 一 2018 年 平均 气温 为 -5.33 C, 
年 均 最 大 .最 小 气温 分 别 为 2.24 9C, 710.74 C, FF 
均 降 水 量 为 383.9 mm, 2012 年 平均 气温 .最 大 气温 
与 最 小 气温 为 最 低 值 。 年 平均 .最 大 、 最 小 气温 变 
化 均 比 较 平 缓 ,降水 量 在 2013 年 之 前 整体 呈 下 降 的 
趋势 ,但 在 2013 年 之 后 呈 逐 年 上 升 的 趋势 (图 2)。 
老虎 沟 大 本 和 音 气 象 站 观测 显示 ,5 月 逐渐 出 现 
日 平均 气温 大 于 0 的 天 数 ,10 月 日 平均 气温 基本 


3.1 
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2 老虎 沟 流 域 气温 、 降 水 量变 化 


Fig.2 Variations of the temperature and precipitation in Laohugou Basin 
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均 小 于 0% ,因此 在 本 研究 区 内 主要 的 消融 季 为 5 一 
9 月 ,7 一 8 月 为 强 消融 期 。 对 于 消融 季 来 说 ,日 平 
均 气 温 为 3.05 Y ,日 最 大 气温 是 12.22 "C ,最 小 气 
温 是 -3.53 % 。 消 融 季 5 一 9 月 月 平均 气温 分 别 
为 -0.50 C 3.55 % .6.00 C.5.53 *C,0.50 ;最 大 气 
温 分 别 为 7.46 C, 10.73 CT. 12.93 C 、12.75 C., 
8.34 C ;最 小 气温 分 别 为 -7.05 C, -2.11 90,0.66 C, 
0% -4.98%。 消 融 季 5 和 -9 月 平均 降水 量 为 39.3 mm, 
82.8 mm 79.6 mm ,60.7 mm,36.3 mm, 占 全 年 总 降水 
量 的 77.5%。7 月 之 前 气温 呈 逐 渐 上 升 的 趋势 且 在 
7 月 达到 最 大 值 ,7 月 之 后 逐渐 减 小 ;降水 量 在 6 月 
之 前 呈 增 加 的 趋势 上 且 在 6 月 达到 最 大 值 ,6 月 之 后 逐 
渐 减 小 (图 2)。 
3.2 径流 量 与 气象 要 素 之 间 的 关系 

为 厘清 气象 要 素 对 冰川 区 径流 的 影响 ,对 径流 
量 与 气温 、 降 水 量 等 指标 的 相关 性 进行 了 分 析 。 图 
3a AEM tir O~3 d 径 流量 的 自 相 关 , 径 流量 的 自 相 关 
系数 在 0.82~0.93 之 间 变 化 且 随 着 时 滞 的 增加 呈 减 
小 的 趋势 。 较 高 的 径流 量 自 相 关系 数 表明 ,在 冰川 
区 冰川 作为 “固体 水 库 ? 对 径流 量 的 响应 起 主导 作 
用 , 即 某 一 天 的 径流 量 与 前 一 天 的 径流 量 高 度 相 
关 。 图 3b~f 为 径流 量 与 气象 要 素 之 间 的 相关 性 ,分 
析 表 明 :径流 量 与 气温 和 水 汽 压 在 时 滞 为 0~3 d 时 


有 较 强 的 相关 性 ,与 气温 的 互相 关系 数 在 0.76~0.86 
之 间 变 化 ,气温 升 高 , 带 来 大 量 的 冰川 融 水 ,气温 是 
控制 冰川 区 径流 的 主要 因素 。 与 水 汽 压 的 互相 关 
系数 在 0.76~0.81 之 间 变 化 ,水 汽 压 是 空气 中 水 汽 所 
产生 的 分 压力 ,间接 表示 大 气 中 水 汽 含量 ,大 气 中 
水 汽 含量 较 大 ,水汽 压 就 较 大 。 水 汽 作为 重要 的 温 
室 气 体 之 一 ,可 以 通过 吸收 太阳 短波 辐射 和 长 波 辐 
射 及 发 射 长 波 辐射 改变 地 气 系统 的 辐射 收 支 ,进而 
影响 地 气 系统 ”。 大 气 中 水 汽 压 增加 ,意味 着 空气 
中 水 汽 含量 增加 ,会 吸收 更 多 的 地 面 长 波 辐射 并 发 
射 向 地 面 , 使 得 地 表 增 温 ,从 而 增加 融 水 径流 , 故 径 
流量 与 水 汽 压 相关 性 较 大 。 与 相对 湿度 的 互相 关 
系数 在 0.42~0.46 之 间 变 化 ,相对 湿度 是 空气 中 的 水 
汽 压 与 相同 温度 下 饱和 水 汽 压 的 百分比 , 即 相 对 湿 
度 也 是 通过 水 汽 压 间接 影响 径流 量 , 因 此 相关 性 较 
水 汽 压 弱 。 径 流量 与 降水 量 相关 性 较 弱 , 与 降水 量 
互相 关系 数 在 0.12~0.27 之 间 变 化 , 即 随 着 降水 的 增 
加 ,径流 量 只 是 微弱 的 增加 ,冰川 上 的 降水 主要 以 
雪 的 形式 发 生 ,反照 率 增加 ,从 而 抑制 冰川 的 消融 ， 
使 得 消融 量 减 少 。 如 果 是 液态 降水 ,其 转化 为 径流 的 
量 小 于 气温 升 高 导致 的 冰川 消融 对 径流 的 补给 
此 相关 性 较 弱 。 与 短波 辐射 无 相关 性 , 随 着 时 灌 的 
增加 ,除了 短波 辐射 ,其 他 要 素 与 径流 量 的 相关 性 


一 一 相关 系数 ---- 当 黑 线 落 在 蓝 色 虚线 内 时 不 相关 
(a) 径流 量 (b) 气温 (c) 降水 量 
1.0 1.0 0.4 
E ES RS 
N in 0.5 ~ 
p p p 
0.0 
3 2 1 0 1 2 3 
ApH /d 时 请 /d 
A KAE (e) 相对 湿度 E 
1.0 0.6 T. (f) 短波 辐射 
六 0.5 M 3 N 
0.0 
3 2 -1 0 1 2 3 
时 滞 /d 


图 3 径流 量 自 相 关 及 与 其 他 气象 要 素 之 间 的 相关 性 


Fig.3 Runoff autocorrelation and correlation with other meteorological elements 
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随 着 时 滞 的 增加 均 有 不 同 程度 的 降低 。 总 体 来 看 ， 
径流 量 与 气温 水汽 压 及 相对 湿度 呈正 相关 , 与 降水 
量 呈 微弱 的 正 相 关 , 这 与 Kang 等 小 在 天 山 及 Li 等 ” 
在 Parlung 4 号 冰川 的 研究 结果 一 致 。 

由 于 在 部 分 年 份 消融 季 径 流 数 据 短缺 ,因此 ， 


与 气象 要 素 建 立 的 回归 方程 所 模拟 的 径流 量 精度 
(图 4)。 经 比较 ,各 年 建立 的 方程 模拟 径流 量 效果 
在 后 2a 较 差 , 多 年 平均 状态 建立 的 方程 模拟 径流 量 
效果 同样 在 后 2a 较 差 ,主要 是 因为 随 着 消融 季 气 温 
升 高 , 冰 下 排水 增加 ,干扰 了 径流 量 与 气象 要 素 之 


本 研究 建立 了 径流 量 与 气象 要 素 之 间 的 多 元 回归 
方程 来 重建 缺失 部 分 的 径流 。 人 参考 相关 研究 ,将 气 
温水 汽 压 作为 自 变 量 得 到 多 元 指数 非 线 性 回归 方 
程 可 以 较 好 地 模拟 冰川 区 径流 3, 故 本 人 研究 将 径流 
作为 因 变 量 ,气温 、 降 水 量 、 水 汽 压 作为 自 变 量 时 所 
得 到 多 元 指数 非 线 性 指数 回归 方程 相关 性 最 高 , 增 
加 降水 量 后 ,模拟 的 效果 较 好 。 本 研究 分 别 建立 了 
2012,2013,2015,2016,2018 年 不 同年 份 消融 季节 
及 5 个 消融 季 平 均 状 态 下 的 多 元 指数 非 线性 指数 回 
归 方 程 ( 表 1) ,由 于 2009 年 仅 有 气温 和 降水 量 的 观 
测 , 故 将 气温 和 降水 量 作为 自 变量 建立 多 元 线性 回 
归 方 程 。 结 果 发 现 ,2009 年 由 多 元 线性 回归 方程 得 
到 的 尼 均 小 于 由 多 元 指数 非 线性 回归 方程 的 尼 。 

由 多 年 平均 径流 量 与 气象 要 素 建立 的 多 元 指 
数 非 线性 回归 方程 模拟 2012.2013 ,2015 、2016 年 和 
2018 年 各 年 消融 季 的 径流 量 ,各 消融 季 的 纳什 效率 
系数 分 别 为 0.68、0.60、0.40、-0.10 及 0.10, 平 均 为 
0.34; 相 对 误差 分 别 为 -1.2%、19.1%、-15.2% 、49.6% 
及 38.2% ,平均 相对 误差 为 18.1%。 相 反 , 由 2012、 
2013、2015、2016 年 及 2018 年 各 年 建立 的 多 元 指数 
非 线 性 回归 方程 模拟 的 各 消融 季 径 流量 的 纳什 效 
率 系 数 分 别 为 0.71、0.94、0.57、0.67 及 0.56, 平 均 为 
0.70 ; 相对 误差 分 别 为 -1.3%、1.5% 、0.5% 、2.2% 
及 -12.4% ,平均 相对 误差 为 -1.9%。 由 此 可 知 ,各 年 
回归 方程 模拟 径流 量 的 精度 高 于 多 年 平均 径流 量 


间 的 关系 (河道 改道 , 测 流 误差 增加 )。 此 外 ,模拟 
峰值 误差 较 大 , 即 模拟 的 峰值 比 实测 的 偏 小 (图 4)， 
主要 原因 是 由 降水 量 、 水 汽 压 造成 的 ,本 研究 利用 
大 本 营 气 象 站 数据 作为 输入 (大 本 营 数 据 时 间 按 序 
列 长 旦 比较 连续 ) ,大 本 营 气 象 站 海拔 为 4180 m, rfi 
冰川 海拔 范围 约 为 4250~5483 m, 相 关 研 究 表 明 老 
虎 沟 流域 降水 量 随 着 海拔 的 增加 而 增加 * ,冰川 的 
降水 量 远 大 于 大 本 营 降 水 量 , 同 时 ,冰川 上 的 水 汽 
压 也 大 于 大 本 营 人 处 的 水 汽 压 , 故 采用 的 降水 量 、 水 
汽 压 输 入 偏 小 ,模拟 的 结果 也 偏 小 ;此 外 , 冰 下 排水 
系统 未 考虑 ,使 得 断面 流量 偏 小 ; 由 于 本 研究 采用 
的 是 统计 关系 来 模拟 计算 日 径流 量 ,不 确定 性 与 系 
统 误差 导致 径流 量 与 实测 有 差异 。 
3.3 径流 量 的 年 际 、 季节、 日 变化 

根据 重建 的 部 分 缺失 径流 ,结合 1959 年 与 1960 
年 的 观测 数据 与 现 有 数据 ,分 析 了 实际 观测 径流 在 
消融 季 的 年 际 、 季 节 和 日 变化 ,尽管 水 文 断面 位 置 
略 有 变动 ,但 基本 维持 在 现 有 断面 附近 ,人 处 于 同一 
海拔 位 置 , 汇 水 面积 变化 不 大 。 此 外 ,冰川 末端 观 
测 的 径流 包括 冰川 消融 区 积累 区 和 裸露 山坡 产生 
的 所 有 径流 。 当 裸露 山坡 面积 在 冰川 区 所 占 比例 
较 小 时 ,冰川 区 径流 与 冰川 融 水 径流 相当 。 所 以 在 
不 同时 间 观 测 的 径流 量 可 以 用 来 分 析 时 间 序 列 的 
变化 。1959、1960、2009、2012、2013、2015、2016 年 
及 2018 年 消融 季 的 平均 径流 量 分 别 为 1.45 ms 


dl 各 消融 季 及 其 平均 状态 下 多 元 回归 方程 的 建立 


Tab. 1 Establishment of multiple regression equation in each ablation seasons and its average state 


年 份 方程 a b c d R 

2009 y=0.501+0.40p - 0.50 - 0.40" 0.51 
2012 y-0.44 e 2er 0.44" 0.24 0.027 1.00 0.71 
2013 y=0.46e 020510090 0.46" 0.20" 0.977 0.003 0.94 
2015 y-1.17e929-971000) 1177 0.237 -0.71 0.01 0.57 
2016 y-0.30 0521000 0.30' 0.15" 2.17" 0.02" 0.67 
2018 y=0.1 50272001» 0.15 0.19” 2o -0.01 0.58 
年 均 y=0.28 00417840020) 0.28 0.24°" 1.797 0.02" 0.89 


UE RON eg ML t 9 Ul sp HEKE so KTR ia 为 方程 的 系数 ;cd 分 别 为 各 个 变量 的 系数 。* \*# x*** 分 别 表示 显著 性 水 平 P<0.05 .P< 


0.01 .P<0.001。 
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图 4 径流 量 观测 值 与 模拟 值 对 比 


Fig.4 Comparison between the measured runoff and simulated runoff 


1.85 m$'*$',2.78 m'*8',1.80 m^*8!, 1.62 m’+s", 
2.30 m*8',224 ms K 1.84 ms1。 产 流 分 别 为 
529 mm,676 mm,1015 mm,657 mm,591 mm,840 mm, 
818 mm,671 mm。1959 一 1960 年 的 平均 径流 量 为 
1.65 ms 2009 一 2018 年 的 平均 径流 量 为 2.10 ms 
可 见 ,21 世 纪 消 融 季 的 平均 径流 量 大 于 20 世 纪 50 
年 代 末 。 为 更 好 理解 径流 量 的 年 际 变化 ,采用 变异 
系数 来 量化 各 消融 季 月 均 径 流量 的 年 际 变化 ,5 一 9 
月 各 月 的 径流 量变 异 系数 分 别 为 0.19、0.27、0.30、 


112.2 mm,269.4 mm .253.4 mm,43.4 mm。 就 消融 下 
各 月 的 产 流 来 说 ,7 月 产 流 最 大 , 占 消 融 季 总 产 流 的 
37.3% ,其 次 是 8 月 , 约 占 35.1%,6 月 约 占 16.1%,5、9 
月 分 别 约 占 5.3%、6.2%。 相 关 人 研究 也 表明 ,对 于 大 
陆 性 冰川 其 径流 量 的 最 大 值 出 现在 7、8 月 ,两 月 的 
径流 量 占 全 年 径流 量 约 为 73%5 ,本 研究 结果 为 
72.4% , HIE ES REA 3X 

径流 量 的 日 变化 特征 :图 6 描述 了 各 消融 季 各 


0.27 及 0.17,6 一 8 月 径流 量变 化 较 大 。 消 融 季 变异 
系数 平均 值 为 0.24 ,表明 整个 消融 季 月 径流 量 的 年 
际 变化 较 小 。 

径流 量 的 季节 变化 :径流 量 从 5 月 开始 增加 ,在 
7 月 达到 最 大 值 ,之 后 开始 减 小 ,8 一 9 月 减 小 的 速率 
较 快 (图 5)。5 一 9 月 月 均 产 流 分 别 为 38.6 mm, 


径流 量 日 变化 特征 。 从 图 中 可 以 看 出 各 消融 季 5 
径流 量 开始 增 加 ,7、8 月 径流 量 达到 最 大 值 ,9 月 
流量 逐渐 减少 直至 消失 。 在 不 同 消融 季 的 开始 
和 结束 , 即 5 月 和 9 月 径流 量 日 变化 都 比较 小 ,径流 
曲线 是 较 平 缓 的 。 在 消融 季 的 开始 ,径流 量 的 日 波 
动 比较 小 ,主要 是 由 于 较 低 的 气温 及 季节 性 积 雪 造 
成 。 此 时 的 融 水 主要 是 积 雪 融 水 且 水 量 比较 少 , 容 
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Fig. 5 Variation of runoff in the ablation seasons 
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图 6 各 消融 季 径 流量 日 变化 


Fig.6 Diurnal variation of runoff in each ablation seasons 


易 被 积 雪 及 冰川 内 部 的 管道 所 拦截 ,此 外 , 融 水 在 。” 长 ,5、6 月 径流 量 达到 最 小 值 与 最 大 值 的 时 间 分 别 
积 雪 层 内 发 生 再 冻结 。 因 此 ,在 消融 季 的 开始 , 径 。 为 06:30 与 15:50、06:20 与 17:00, 可 见 6 月 径流 量 到 
流量 变化 比较 小 (5、6 月 的 标准 差分 别 为 0.20、 ， 达 末 端 时 间 有 所 延长 。 

0.46) 。 径 流 由 冰川 表面 到 达 末 端 经 历 的 时 间 较 在 强 消 融 期 径流 量 的 日 变化 逐渐 变 大 (7、8 月 
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的 标准 差分 别 为 1.00.1.43) ,尤其 在 7.8 月 径流 量 上 
升 、 下 降 均 比较 快 。 主 要 原因 是 气温 升 高 .太阳 辐 
射 增强 使 得 融 水 增加 (冰川 融 水 为 主 )。 冰 川 融 水 
受气 温 影响 比较 大 ,由 于 昼夜 温差 较 大 所 以 径流 量 
日 变化 也 比较 大 。 另 外 ,降雪 量 的 减少 与 冰川 表面 
水 系 的 发 展 降低 了 冰川 对 融 水 的 拦截 作用 ,使 得 融 
水 的 日 变化 更 强 。7、8 月 径流 量 达到 最 小 值 与 最 大 
值 的 时 间 分 别 为 07:20 与 15:50、07:20 与 15:20。 相 
比 6 月 答 流 到 达 末 端 时 间 有 所 提前 。 

至 消融 期 末 , 即 9 月 , 随 着 气温 的 降低 及 降雪 的 
RÆ ,冰川 融 水 减少 ,径流 量 的 日 变化 相 比 强 消 融 
期 大 幅 减 小 (标准 差 仅 为 0.28)。 虽 然 此 时 冰川 上 降 
雪 量 有 所 增加 ,但 积 雪 会 在 很 短 的 时 间 内 消融 ,万 
其 当天 气 状况 比较 好 (晴天 ) 的 时 候 , 径 流 到 达 末 端 
的 时 间 会 比 强 消融 期 有 所 提前 ,主要 是 由 于 较 发 达 
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的 冰 层 表面 水 系 。 此 时 ,径流 量 达 到 最 小 值 与 最 大 
值 的 时 间 为 07:00 与 15:00。 
3.4 径流 的 时 滞 效 应 

冰川 径流 的 时 滞 指 日 内 气温 达到 最 大 值 与 径 
流量 达到 最 大 值 之 间 的 时 间 间 距 。 图 7 为 消融 季 各 
月 气温 及 各 月 径流 量 的 日 变化 ,并 标注 了 最 大 气温 
(7 与 最 大 径流 量 (O.) 出 现 的 时 间 。 老 虎 沟 流域 
径流 具有 显著 的 日 变化 过 程 ,流量 曲线 呈现 * 谷 - 
峰 ” 的 日 变化 特征 。 在 消融 期 初 ,5 月 径流 量 达 到 
最 大 值 (0 ) 的 时 间 述 于 气温 达到 最 大 值 (7 ) 的 
时 间 , 尤 其 在 6 月 ,0 出 现时 间 (17:00) 远 迟 于 7 
(14:00) 且 差距 最 大 ,时 浦 为 3h, 主 要 原因 是 6 月 降 
水 量 较 大 且 降 雪 量 比重 大 ,产生 的 融雪 径流 易 受 到 
积 雪 的 拦截 ,另外 ,消融 期 初 冰川 表面 的 水 系 不 发 
达 , 所 以 融雪 径流 到 达 末 疝 的 时 间 被 延长 0, 与 
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图 7 消融 季 各 月 径流 量 与 气温 日 变化 


Fig. 7 Monthly runoff and temperature diurnal variation in the ablation seasons 
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7 之 间 的 时 间 间 距 最 大 。 

在 强 消融 期 ,降雪 量 减 少 , 降 水 多 以 液态 的 形 
式 发 生 , 雪 的 阻碍 作用 降 至 最 低 。 随 着 气温 的 升 
高 ,冰川 融 水 的 增加 ,使 得 冰川 表面 水 系 逐 渐 发 
展 。 故 虽然 7 出 现时 间 (15:00) 相 比 消融 期 初 有 所 
推迟 ,但 0, 出 现时 间 有 所 提前 ,7 与 0 出 现时 间 
较 消 融 期 初 接近 ,两 者 基本 一 致 。 此 时 ,径流 量 与 
气温 的 变化 曲线 相似 ,这 与 i 等 3 的 研究 结论 相 
似 , 即 在 强 消融 期 两 者 出 现 的 时 间 较 为 一 致 。 到 达 
消融 期 末 , 即 9 月 温度 逐渐 降低 ,降水 多 以 雪 的 形式 
发 生 , 冰 川 径流 逐渐 减少 ,加 上 雪 的 拦截 作用 ,虽然 
Ta 出现 时间 (14:00) 相 比 强 消融 期 有 所 提前 ,但 Qna 
出 现时 间 (15:00) 与 强 消融 期 基本 相同 ,7 与 0% 出 
现 的 时 间 间 距 较 强 消融 期 有 所 增加 , 介 于 消融 期 初 
和 强 消融 期 之 间 。 


4 讨论 


冰川 末 问 观测 的 径流 包括 冰川 消融 区 、 积 累 区 
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38 , Uy 0.46; 与 降水 量 的 正 相 关 性 最 弱 , DU 0.27 
(图 3)。 故 冰川 融 水 径流 受气 温 变 化 更 显著 ,而 降 
水 量 影响 相对 较 小 ,这 与 刘 泵 更 等 2 在 天 山 乌 鲁 木 
齐 河 源 1 号 冰川 的 研究 结论 一 致 。20 世 纪 0 年 代 
未 消融 季 平 均 径 流量 为 1.65 ms 而 21 世 纪 消 融 
季 平 均 径 流量 为 2.10 ms ,主要 受气 温 升 高 影响 ， 
对 应 气温 20 世 纪 50 年 代 来 消融 季 平 均 气 温 为 
1.35% ,而 21 世纪 消融 季 平 均 气 温 为 2.10 % ,消融 
季 气 温 上 升 0.75 Co 

消融 季 强 消融 期 (7 一 8 月 ) 产 流 分别 占 整个 消 
MÆ HI 37.396 , 35.196 , 主要 是 因 强 消融 期 气温 为 
5.8 % ,水 汽 压 平均 为 0.55 kPa, 相 对 湿度 为 64.4% , 
降水 量 占 消融 季 的 50.8%, 较 高 的 气温 使 得 冰川 产 
生 大 量 融 水 ,此 外 ,水 汽 压 与 相对 湿度 较 大 ,降水 量 
较 大 也 使 产 流 有 所 增加 。6 月 产 流 占 消融 季 的 
16.1% ,径流 主要 是 融雪 径流 ,6 月 平均 气温 为 
3.6 % ,降水 量 占 消融 季 的 27.7% ,根据 陈 仁 升 等 " 
的 西北 高 寒山 区 两 雪 分 离 温度 为 3.5~5.5% ,6 月 降 


和 裸露 山坡 产生 的 所 有 径流 。 当 裸露 山坡 面积 在 
冰川 区 所 占 比 例 较 小 时 ,冰川 区 径流 与 冰川 融 水 径 
流 相 当 下。 本 研究 断面 设 定 在 距 末 端 约 1 km Ah, Wt 
面 观 测 到 的 流量 包含 裸露 山区 的 产 流 , 但 由 于 面积 
较 小 ,所 以 观测 的 径流 主要 为 冰川 融 水 径流 。 在 消 
融 季 ,冰川 融 水 径流 的 变化 主要 受气 象 要 素 影响 。 
径流 量 与 气温 为 正 相 关 , 且 相关 性 最 高 (0.86); 与 水 
汽 压 正 相关 ,相关 性 为 0.81; 与 相对 湿度 正 相 关 性 较 


雪 量 占 降水 总 量 的 范围 为 64.4%~82.9%。 所 以 虽然 
6 月 降水 量 较 大 ,加 上 主要 以 降雪 的 形式 发 生 (图 
8) ,但 气温 也 较 强 消融 期 低 , 故 产 流 较 7、8 月 次 之 。 
在 消融 期 初 、 末 (5、9 月 ), 产 流 分 别 占 消融 季 的 
5.3% 6.2% , 主要 是 因 5 .9 月 气温 较 低 ,平均 气温 分 
别 为 -0.18 % .0.33 %。 虽 然 5 月 的 降水 量 (39.3 mm) 
大 于 9 月 (36.3 mm) (主要 以 降雪 形式 发 生 ), 但 5 月 
的 气温 低 于 9 月 ,所 以 产 流 略 低 于 9 月 。 整 体 来 看 5 
月 与 9 月 由 于 较 少 的 降水 量 及 较 低 的 气温 ,导致 产 


100 BS 降水 量 DD 降雪 量 一 4 一 气温 8 
E 4 
SB ZH 
E 27 
= 
B 
0 
-2 
7 8 9 
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图 8 消融 季 降 水 量 降雪 量 与 气温 的 分 布 


Fig. 8 Distribution of precipitation/snowfall and temperature in the ablation seasons 
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大 值 的 时 间 远 迟 于 气温 达到 最 大 值 的 时 间 ,尤其 在 
6 月 ,主要 因 消 融 期 初 降雪 量 较 大 ,5 月 均 为 降雪 ,6 
月 降雪 量 为 64.4%( 雨 雪 分 离 温 度 为 3.5 C)( 图 8)，,5 
月 由 于 气温 较 低 ,融雪 径流 较 少 ,到 了 6 月 随 着 气温 
的 升 高 , 积 雪 开始 融化 ,但 融雪 径流 容易 受到 积 雪 
的 拦截 ,再 加 上 冰川 表面 水 系 不 发 达 , 所 以 融雪 径 
流 到 达 末 端的 时 间 会 被 延长 ,与 气温 达到 最 大 值 的 
时 间 间 距 最 大 。 相 反 , 在 7.8 月 降雪 量 减 少 ,降雨 量 
增加 ,气温 增高 ,冰川 表面 水 系 发 达 , 所 以 径流 到 达 
末端 时 间 会 缩短 ,与 气温 达到 最 大 值 时 间 间 距 最 
小 ,尤其 在 8 月 最 明显 。 所 以 ,冰川 区 径流 除了 受气 
温水 汽 压 、 相 对 湿度 .降水 量 等 因素 的 影响 ,还 会 
受 积 雪 的 影响 ,对 冰川 径流 的 研究 应 综合 考虑 上 述 
因素 的 影响 。 


5 结论 


(1) 老虎 沟 流 域 气温 在 7 月 之 前 呈 增 加 趋势 ， 
之 后 为 减 小 趋势 ,降水 量 在 6 月 之 前 为 增加 趋势 ,之 
后 减少 ,近年 来 气温 呈 上 升 趋势 ,降水 量 亦 有 所 增 
加 。 对 于 径流 量 ,其 自 相关 系数 在 0.82~0.93 之 间 变 
化 ,有 较 强 的 自 相 关系 数 , 与 气温 相关 性 最 高 
(0.86) ,与 水 汽 压 相 关 性 为 0.81, 与 相对 湿度 为 
0.46, 与 降水 量 相关 性 最 弱 仅 为 0.27, 可见 径流 量 受 
气温 影响 最 大 。 由 各 年 径流 量 与 气温 水汽 压 及 降 
水 量 建 立 的 多 元 指数 非 线 性 回归 方程 可 以 较 好 地 
模拟 各 年 的 径流 量 , 多 年 平均 纳什 效率 系数 为 0.70， 
平均 相对 误差 为 -1.9%。 

(2) 21 世纪 平均 径流 量 为 2.10 m^ s", $E 20 tlt 
2 50 FARR ÉY 1.65 ms 有 所 增加 ,主要 因 消 融 季 
气温 增加 0.75 CC 所 致 。 强 消融 期 径流 量 年 际 变化 
较 大 ,消融 期 初 和 末年 际 变化 较 小 。 在 消融 季 , 径 
流量 从 5 月 开始 增加 ,在 7 月 达到 最 大 值 , 之 后 开始 
减 小 ,5 一 9 月 月 产 流 分 别 占 整个 消融 季 的 5.3% 、 
16.1% 、37.3% .35.1% .6.2%。 对 于 消融 季 各 月 径流 
的 日 变化 ,消融 期 初 和 末 ,径流 的 日 变化 比较 小 , 径 
流量 到 达 最 小 值 与 最 大 值 的 时 间 间 距 较 大 ,但 在 强 
消融 期 ,径流 的 日 变化 较 大 ,径流 量 达 到 最 大 值 与 
最 小 值 的 时 间 间 距 较 小 ,主要 受气 温 \ 降 水 / 雪 量 水 
汽 压 、 相 对 湿度 、 冰 川 表面 水 系 等 因素 的 影响 。 

(3) 径流 的 时 湿 效 应 ,老虎 沟 流 域 消融 季 各 月 


径流 量 呈 现 “ 谷 - 峰 ” 的 日 变化 特征 ,消融 期 初 气 温 
达到 最 大 值 与 径流 量 达 到 最 大 值 的 时 间 间 距 较 大 ， 
尤其 在 6 月 两 者 差 3h( 气 温 达 到 最 大 值 时 间 为 14:00， 
径流 量 达 到 最 大 值 时 间 为 17:00) , 主要 受降 雪 量 A 
温 、 冰 川 表 层 水 系 等 的 影响 。 在 强 消融 期 两 者 的 时 
间 间 距 最 小 ,径流 量 与 气温 的 日 变化 曲线 相似 , 主 
要 因 气 温 升 高 ,降水 量 增加 ,冰川 表面 水 系 发 达 。 
消融 期 末 两 者 关系 介 于 消融 期 初 与 强 消 融 期 之 间 ， 
主要 因 降 雪 比 重 增加 。 
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Abstract: Based on the runoff data and meteorological data of Laohugou watershed in the western of Qilian 
Mountains, northwest China, the research analyzed the relationship between runoff and meteorological variables 
in glacier area, and established a multivariate exponential nonlinear regression to reconstructe the runoff. In 
addition, the characteristics of interannual, seasonal, diurnal variations of runoff in the glacial region were 
analyzed. The results showed that: (1) The correlation between runoff and temperature is the highest (0.86), 
followed by water vapor pressure (0.81), relative humidity (0.46), and precipitation (0.27). Runoff is most 
affected by temperature. (2) The mean daily runoff is 2.10 m°» s™' in this century, which is higher than 1.65 m?^:s'' 
in the late 1950s, the main reason is that the temperature increased by 0.75 °C in the ablation season. The 
interannual variation of runoff is large in the strong ablation period, and the interannual variation is small in the 
beginning (May-June) and end (September) of the ablation period. The proportion of runoff generation from May 
to September was 5.3%, 16.1%, 37.3%, 35.1% and 6.2%, respectively. (3) Multivariate exponential nonlinear 
regression equation can better simulate the daily runoff (the mean Nash efficiency coefficient is 0.70). After the 
supplement of the missing runoff, the diurnal variation of runoff is small in the beginning and end of the ablation 
period, but the diurnal variation of runoff is large in the strong ablation period. For the time-lag effect of runoff, 
monthly runoff in Laohugou watershed showed a characteristic of diurnal variation of “valley-peak” during the 
ablation period. The time interval between the maximum temperature and the maximum runoff was long in the 
beginning and end of the ablation period, while the time interval was short between the maximum temperature 
and the maximum runoff during the strong ablation period, and the maximum difference between the maximum 
temperature and the maximum runoff was 3 hours in June. 
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